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ABSTRACT 

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for 
a new sample of simulated field galaxies using fully cosmological, ACDM, high resolution SPH + 
N-Body simulations. These simulations include an explicit treatment of metal line cooling, dust and 
self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos 
have masses ranging from a few times 10 8 to nearly 10 12 M Q . At z=0 our simulated galaxies have a 
baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10 3 - 4.5 
x 10 10 Mq we find extremely good agreement between the SHM ratio in simulati ons and the present - 
day predictions from the statistical Abundance Matching Technique presented in lMoster et al.l (|2012ft . 
This improvement over past simulations is due to a number systematic factors, each decreasing the 
SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold 
gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and 
photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 
30% overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over 
cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes 
can underestimate the contribution of old stellar populations to the total stellar mass, leading to 
stellar mass errors of up to 50% for individual galaxies. These results highlight the importance of 
using proper techniques to compare simulations with observations and reduce the perceived tension 
between the star formation efficiency in galaxy formation models and in real galaxies. 
Subject headings: galaxies: evolution — galaxies: formation — methods: N-Body simulations 
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19801 iBlumenthal et all 119841: iDekel k Silfl 119861 : 
White k Frenkl I 1991f ). many galaxy properties are 



expected to correlate with the mass of the galaxy's host 
halo. In particular, the stellar-to-halo mass relation 
(SHM), defined as the ratio of the stellar mass (M stor ) 
within a halo of total mass M-halo within a given 
over-density (< p > / p cr it = 200 in this work) is a 
robust estimator of the efficiency of gas cooling and 
star for mation (SF) processes over a wide range of hal o 
masses (|Somerville k Primac afl99l lB owcr "etaI|[20T0h. 



Both observational (|Heavens et al 
2007) and theoretical work 



~04[ iZhengetai] 
(jBower et al.l 120111 ) suggest 



that the SF efficiency peaks at the scale of L* galaxies 
and declines at smaller and larger masses. On the 
low mass end, this decline is likely because SF is 
suppressed by gas heating from the UV cosmic field 
(jQuinn et al.l I1996I; iGnedinl 120001; lOkamoto etaLl [20081: 
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iNickerson et~aT1l2011[ ) and supernova (SN) heating with 
gas removal. At larger masses, energy feedback from 
super-massive black holes (SMBHs) is thought to be the 
domina nt process responsible for lowering the SF effi- 
cienc y dBower et al.ll2006t lCrotoiill2009t McCarthy et all 
120111 1 Johansson et alJl2012h . 

Recently. iMoster et all (120121 hereafter M12) and 
other groups flVale k Ostrikerll200i iConroy et alj|2006t 



M 
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Mandelbaum et al.l 120061: iMore et all 120091: IGuo et al I 



2010t iTruiillo-Gomez et all 120111: IBehroozi et al.l 12012ft 

used the Abu ndance Matching Technique (AMT) and 
its variations (Ya ng et al.l [2012) to derive a SHM re- 
lation of real galaxies. In its simplest form AMT as- 
sumes a monotonic relation between the stellar mass 
function of (real) galaxies and the underlying halo mass 
function. This relation is constrained by matching the 
observed galaxy stellar mass function to the ACDM 
halo mass function from N-body simulations. Simi- 
lar works ([Guo et al.ll2oTot iLeauthaud et al.ll20lU I2012t l 
have included constraints from lensing and used slightly 
different underlying cosmologies. This approach has 
also been used to c onstrain the scatter in the SHM 
([Reddick et al.l [20121 ) by c omparing the predicted spa- 
tial c l ustering of DM halos dSheth et all2001t iReed et all 
120071 Ivan Daalen et al.l 120121 ) with the observed abun- 
dances and clustering properties of galaxy population s 
(IBlain et al1l200"i IConrov et all 120011 IReid et all[2010l) . 
Additionally. IBehroozi et all (|2012f) ~discuss the implica- 
tions of the upturn in the faint-end slope of the stellar 
mass function on the SHM relationship. 

Several works have highlighted how uncertainties in 
the derived SHM relation depend on a number of fac- 
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ters, some of them poorly known. For example the 
stellar masses of real galaxies are inferred from opti- 
cal and near-I R photometric measur e ments and/or re- 
solve d spectra (jBell k de Jong) 120011 : iKauffmann et al.l 
2003). This approach carries substantial unc e rtain- 
ties and possib l e degenerac i es dBell k de Jong] 1 2001 



Bell et al.l [2001 IPforr et al.l 1201% IHuang et al.l 12012 
Behroozi et al.l I2010l) as the observed spectral energy 



distribution of a galaxy is a function of many phys 
ical processes (e.g., stellar evolution, SFH and metal 
enrichment history, and wavel ength- dependent dust at- 
tenuation, iPanter et al.l I2004D . Furthermore, as sur- 
veys often measure individual galaxy magnitudes within 
an aperture based on a surface brightness cutoff, the 
mass of the stellar component co uld be systematically 
underestimated by at least 20% (|Graham et al.l 120051 : 
Shimasak u et al.l [20011 ) if part of it is old (hence faint) 
and/or low surface brightness. Finally, the number den- 
sity of galaxies will be affected by incomplete ness at the 
faint end of the galaxy luminosity function (IDalcanton 
19981: ISawala et al.ll2011trGeller et al.H20T2 : Santin i et al 



20121) . Separate from worries over the stellar mass 
determinations are worries about the halo mass func- 
tion. While the halo mass function ob tained in DM-only 
simulations is robu stly constrained (IReed et al.l 120031 : 
Sprin geTet al.|[2005l ). it has recently ()Sawala et al.ll2012T) 
been shown that halo masses in DM-only simulations 
exceed those obtained in simulations including baryon 
physics by up to 30%, introducing another systematic 
bias in the AMT, as a galaxy of a given stellar mass 
is matched with a too massive halo, pushing the stellar- 
halo mass ratio down. Taken together, the above caveats 
suggest that a better understanding of the connection 
between galaxy masses and the underlying halo masses 
could in principle be gained by using realistic simulations 
of galaxy formation that directly include baryon physics, 
SF and SN feedback. 

While substantial progress has been made in 
creatin g galaxies from cosmological initial condi- 
tions (IScannapieco et all 120101: [ McCarthy et al.l 120121: 
Sales et all 2012; IStinsorT et al.l 120121 : Uohansson et al.l 



20121) . several recent studies have pointed out a large 



discrepancy between the SHM relation estimated for 
real galaxies and the one obtained in several numerical 
simulations of galaxy formation (jSawala e t al.l 120111 : 
IGuo et all I2010D . Simulations have repeatedly shown 
that a lack of realistic SN feedback leads to overesti- 
mating sta r formation a s part o f the general over c ooling 
problem (lAbadi et aD 120031: iGovernato et al.l 120071 : 



iPiontek k Steinmetd 1201 It iKeres et al.l 1201 1| ). Most 
simulations that overproduce stars f orm galaxies tha t 
have large spheroidal components (|Eke et al.l 120011) . 
Increment al improvements based on mo r e realistic SN 
feedb ack (|Thacker k Couchmanl l2000t iStinson et al.l 
120061 ) led to si mulations that formed galaxies with 
exten ded disks (jGovernato et al.l 120091 : iBrooks et al.1 
120 111 ), but still substantially overproduced stars. Only 
recently, a new generation of high resolution simulations 
demonstrated the impact of feedback at lowering SF 
efficie ncy such that SF occurs only at high gas densi- 
ties (|Ceverino k Klvpinl 120091: IGovernato et al.l 120101: 
Guedes et al.ll2011l: IGovernato et aL| [2012: Zo lotov et al.l 
2012tlBrook et al.ll2012t ). As SF is more efficient in dense 
gas clouds, feedback from these high density regions 



generate outflows that simultaneously improve on several 
long standing proble ms namely the subs t ructure over- 
abundance problem (jMoore et al.l Il998t iKlypin et al.l 
1999; Benson 201C|), reducing the B/D ratio in small 
galaxies by removal of low angular momentum baryons 
dBinnev et al.l 120011: IGovernato et al.ll2010t iBrook et all 
1201 If ) and forming DM cores by transferring en- 
ergy fr om baryons to the DM (Mashc henko et al.1 



20061: iPasetto et all 
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I2012t IGovernato et al.l 12012ft . Forming 



stars in dense gas regions is a crucial step, as observa- 
tions strongly support that the spa tially resolved SF 
is linked to the loca l H 2 fraction (iBigiel et al.l 120081 : 
iKrumholz et~aT1 120091: iGenzel et all |2012|), which only 
becomes significant at the density of star forming 
regions, > 10-100 amu/cm 3 . 

The relationship between gas density and H2 abun- 
dance can now be naturally implemented in sim- 



ulations, and the creation and destruction of H9 
can be followed con sistently ()Gnedin et al.l 12009c 
iChristensen et ah! 12012). allowing much more realistic 
simulations to be run, and establishes a physically moti- 
vated connection between SF and high density (shielded) 
gas. This is indeed one of the major steps forward in 
the work presented here. While feedback from SNe re- 
mains still poorly understood, its effects are being ob- 
served over a large range of reds hifts and galaxy masses 
(Martin 1999: 1 Wang et al.ll2010ft . It is therefore impor- 
tant to evaluate a new set of high resolution simulations 
to test if outflows can form galaxies with realistic obser- 
vational properties that also reside on the SHM relation. 

Relatively less attention has been given to comparing 
results from simulations with observational estimates of 
the SHM relation in a consistent way. While some recent 
work s reporte d ( S awal a et al.l 120111: lAvila-Reese et al.l 
120111: IPiontek k Steinmetzl 120111: iLeitnerl I2012D an ex- 
cess of stars formed in simulations, they compared the 
galaxy stellar and total halo masses directly measured 
from simulations while those quantities are actually in- 
ferred from the light distribution of real galaxies. A con- 
sistent ap proach has been already tried with promising 
results in lOh et al.l (|2011f ). where two simulated dwarfs 
were compared with a set of galaxies from the THINGS 
survey, finding that the simulations have the same central 
baryon and DM distribution as in the observational sam- 
ple. In that work, stellar and halo masses were obtained 
using photometric and kinematic data for both the sim- 
ulat ed and th e real sample, finding excellent agreement 
(sec IQh et all (I20T1 figure 5). 

In this paper, we present a consistent comparison be- 
tween the SHM estimated in Moster et al. (2012) and 
a set of high resolution simulations spanning 5 orders of 
magnitude in stellar masses. The evolution of galaxies 
is simulated at high resolution using Smoothed Particle 
Hydrodynamics (SPH) in a cosmological context. We fo- 
cus on how galaxies populate dark matter halos ranging 
from 10 8 — 1O 12 M in a field environment and comparing 
results with the estimates of the SHM from M12. The 
four most massive, MW-like galaxies have simil ar reso- 
lution to the "Eris" galaxy (jGuedes et al.l l20liT ). while 
the smaller galaxies have even better force resolution ( 
down to 65 pc and star particles as small as 450 M©). 
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Fig. 1. — Baryonic fraction with respect to the cosmic ratio, 
for simulated field galaxies as a function of stellar mass, measured 
at z=0. Circles are the "direct from simulation" results, includ- 
ing all gas and stars within R.200- Triangles are the "observable" 
baryon fractions, includind all stars and all the 'observable' cold 
gas (defined as 1.4 X (HI+H2), within R,20o)- The empty symbols 
are galaxies with no observable gas (cold gas mass < 100 Mq). 
Galaxies below 10 8 Mq lose a significant fraction of baryons due 
to heating from the cosmic UV background and SN feedback. 

All simulations include the effects of metal line cooling 
and H2 dependent star formation. 

Our simulations include cooling, star formation, a cos- 
mic UV background and form galaxies with structural 
properties comparable to the real ones. The dataset uti- 
lized is described in detail in Section 2 and has been ana- 
lyzed in other papers sho wing that the galaxies fol low the 
Kennicutt-Schmidt law ([Christensen et aLl 12012ft. have 
cored DM profiles similar to the observed ones et al.l 
I201H iGoyernato et al.l 12012ft and realistic satellit e pop - 
ulations (jZolotov et al.l l2012t iBrooks fc Zolotovl 12012ft . 
Without any further fine tuning, in this work we com- 
pare the same simulations to the SHM relation obtained 
in Ml 2, using an analysis technique comparable to the 
one used in the original paper to estimate their stellar 
and halo masses. We show that simulations form realis- 
tic galactic systems that also match the z=0 SHM of real 
galaxies over five orders of magnitude in stellar mass. 

The paper is organized as follows: in §2 we describe 
the details of our N-body simulations. In §3 we compare 
results with the SHM predicted in M12. The results are 
discussed in §4. 

2. THE SIMULATIONS 

The simulations used in this wor k were run with the 
N-Body + SPH cod e GASOLINE (|Wadslev et al.ll200i 
iStinson et al J 12006) in a fully cosmological ACDM con- 
text: ft = 0.26, A=0.74, h = 0.73, cr 8 =0.77, n=0.96. 
The galaxy sample was selected from two uniform DM- 
only simulations of 25 and 50 Mpc per side. From these 
volumes a few field-like regions were selected and then 
resimulated at higher resolution u sing the 'zoomed-in ' 
volume renormaliza tion technique ()Katz fe White! IT9931 ; 
iPontzen et al.l 12008ft . This technique allows for signif- 
icantly higher resolution while faithfully capturing the 
effect of large scale to rques that deliver angular m omen- 
tum to galaxy halos (jBarnes fc Efstathioul 11987ft . With 
this approach, the total high resolution sample contains 
eighteen field galaxies, each halo resolved by 5xl0 4 to a 



few 10 6 DM particles within TL„i r , defined as the radius 
at which the average halo density = 200 x p cr n . 

The force spline softening ranges between 64 and 
170 pc in the high resolution regions of each volume 
and it is kept fixed in physical coordinates at z < 10. 
Star particles are formed with a mass of 400-8000 Mq. 
The halo mass range covered by the simulations spans 
nearly four orders of magnitude, from a few times 10 s 
to 8 x 10 11 M Q (peak velocities V pea k = 10 to 200 
km/sec), and stellar masses M sta r from 10 4 to a few 
10 10 Mq. As other works have highlighted the impor- 
tance of having a r epresentative samp l e before drawing 
general conclusions (IBrooks et allfeOllHSales et al.ll2012t 
iMcCarthv et al.l 12012ft . the halos in our sample span a 
represent ative range of ha lo spin values and accretion 
histories (jGeha et al.1 12006). Galaxies and their parent 
halos were first identified using AHE0 (jGill et al.l 120041: 
iKnollmann fc Knebell2009ft . The total halo mass (includ- 
ing DM, gas and stars) is defined at a radius R V i r , defined 
as the radius at which the average halo density = 200 x 
Pcrit, consistent with M12. No sub-halos have been in- 
cluded in our sample, although the most massive galax- 
ies ha ve a realistic population of satellites ()Zolotov et al.l 
12012ft . More details on this dataset a nd the properties 
of the sate llite population are given in iGovernato et al.1 
(l20Tl and IBrooks fc Zolotovl ll20Tl . 

2.1. H2 fraction, Star Formation and SN Feedback 

In a significant improvement, t his new set of sim ula- 
tions include metal line cooling (jShen et al.l 12010ft and 
a dust dependent description of H2 creation and de- 
struction b y Lyman- Werner radiation and shielding of 
HI and H 2 (|Gnedin et all 120091: IChristensen et all 120121 
hereafter CH12). As in CH12, the star formation rate 
(SFR) in our simulations is set by the local gas density 
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Fig. 2. — The stellar mass ratio between the galaxies simulated 
with the old 'low density SF threshold' and the new sample. In the 
new sample SF is regulated by the local abundance of molecular 
hydrogen, resulting in feedback significantly lowering the total SF 
efficiency. All quantities as measured directly from the simulations. 



8 A MIGA'S Halo Finder, available for 
http: / /popia. ft. uam.es/AHF/Download. html 
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ciency parameter, c* = 0.1, gives the correct normaliza- 
tion of the Kennicutt- Schmidt relation (the SF efficiency 
for each star forming region is much lower than the im- 
plied 10%, as only a few star particles are formed before 
gas is disrupted by SN winds) . With the inclusi on of 
the H2 fraction term (see also iKuhlen et all 1201 If ) . the 
efficiency of SF drops to zero in warm gas with T > 
3,000 K. The simulations include a scheme for turbu- 
lent mixi ng that redistribu tes heavy elements among gas 
particles ()Shen et al.ll2010D . With this approach, the lo- 
cal SF efficiency is linked to the local H2 abundance, as 
regulated by the gas metallicity and the radiation field 
from young stars, without having to resort to simplified 
approaches based on a fixed local gas dens ity threshold 
(IGovernato et al.ll20Tot Ikuhlen et al.ll20Tl . The simu- 
lations assumed a Kroupa IMF and relative yields, but 
observable quantities have been co nverted to a Chabri er 
IMF, for a direct comparison with iMoster et all (|2012D . 

As in previous works using t h e "bl a stwave" SN feed- 
back approach (jStinson et al.1 120061 iGovernato et al.1 
2012), mass, thermal energy, and metals are deposited 
into nearby gas when massive stars evolve into SNe. 
The amount of energy deposited amongst those neigh- 
bors is 10 51 ergs per SN event. Gas cooling is then 
turned off until the end of the momentum-conserving 
phase of the SN blastwave which is set by the local 
gas density and temperature and by the total amount 
of energy injected, typically a few million years. Equi- 
librium rates are computed from the phot oionization 
code C loudy (jFerland et al.lll998t ). following [SheneFal] 
(2010). A spatially uniform, time evolving, cosmic UV 
background turns on at z — 9 and modifies the ioniza- 
tion and excitation state of t he ga s, following an updated 
model of lHaardt fc Madaul (|1996h . 

This feedback model differs compared to other 
"sub-grid" schemes (e. g., ISpringel fc Hernquistl 120031 : 
IScannapieco et aLll2012j ) in that it keeps gas hydrody- 
namically coupled while in galactic outflows. The ef- 
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Fig. 3. — The cold gas mass as a function of stellar mass. Sim- 
ulations vs. SHIELD and ALFALFA data. The HI mass of each 
galaxy in the simulated sample is plotted vs the SDSS r-band mag- 
nitude and compared to two samples from nearby surveys. Red 
solid dots: simulations. Diamonds: ALFALFA survey. Asterisks: 
SHIELD survey. While feedback removes a large fraction of the 
primordial baryons, the simulated galaxies have a high gas/stellar 
mass ratio, comparable to the observed samples. Most of the 
cold gas resides within a few disk scale lengths from the simulated 
galaxy centers. 



ficient deposition of SN energy into the ISM, and the 
modeling of recurring SN by the Sedov solution, should 
be interpreted as a scheme to model the effect of energy 
deposited in the local ISM by all processes related to 
young stars, including UV radiation f rom massive stars 
(|Hopkins _et alJl201ltlWise et al.l[2012h . TheSFHsofthe 
galaxies in our simulated sample are bursty, especially 
those res iding in halo s smal ler t han 10 10 M^. As dis 



cussed in IBrook et al.l ([201 ID and lPontzen fc Governatol 
(2012), a bursty SFH is necessary to remove low angu- 
lar momentum baryons and create the fast outflows able 
to transfer energy fr om baryons to the DM component 
and create DM cores (jGovernato et al.ll2010f) . These pro- 
cesses lead to realistic dwarf galaxies with slowly rising 
rotation curves and typical central surface brightnesses 
21 < fiB.o < 23.5 (|Oh et alj|20ll[ ). Outflows and the 
cosmic UV field progressively suppress star formation in 
halos with total mass smaller than a few 10 10 M©. Test 
runs verify that the effect of energy feedback on sup- 
pressing SF is much l arger than that of having a low H2 
fraction-SF efficiency (Christcnscn et al. 2012). 

2.2. Star formation efficiency as a function of galaxy 
halo mass 

As a result of SF feedback processes, the SF efficiency 
is greatly reduced over the whole mass range of our sim- 
ulated sample. The smallest galaxies in our sample turn 
only ^0.01% of their primordial baryons into stars. The 
more massive galaxies in our sample turn ^30% of their 
primordial baryon content into stars, but we demonstrate 
in the next section that using stellar masses based on 
photometry reduces this efficiency to an apparent 10%. 
Feedback expels about 70% of the gas to outside of R V i r 
in dwarf galaxies with v c ~ 40-50 km s _1 . Larger galaxies 
retain a larger fraction of their original baryons, while the 
smallest galaxies lose an even larger fraction of baryons 
due to gas heating by the cosmic UV background, which 
further reduces their SF (see Figure 1). 

To evaluate the effect of the adopted SF model on 
the resulting SF efficiency of our galaxy sample, galax- 
ies were re-simulated using the SF approach used in our 
older works. These reference runs adopt a lower density 
thresh old for SF, 0.1 amu/cm 3 . As discussed in several 
works (|Governato et al.ll2010t IGuedes et al.ll2011l) . a low 
density threshold makes SF more diffuse and locally less 
efficient. In this scenario, typical of low resolution simu- 
lations where star forming regions cannot be resolved, the 
amount of SN energy per unit mass delivered to gas par- 
ticles is effectively lowered, making feedback much less 
effective at suppressing SF. While the gala xies in the low 
thres hold sample have realistic disk sizes (|Brooks et al.l 
1201 If ) and the total amount of energy released into the 
gas is actually a few times larger, they contain many 
times more stars (see Figure 2) and overproduce stars 
compared to the SHM relation. The comparison between 
the amount of stars formed in the old and the new runs 
(Figure 2) demonstrates that the large decrease in the 
SF efficiency in the new simulations (as much as a factor 
of 15) is due to the improved i mplementation of SF i n 
dense/H 2 rich regions (see also IGovernato et a.1.1 12010h . 
This lower SF efficiency goes a long way toward recon- 
ciling simulations of galaxy formation with current esti- 
mates of the SHM relation. 



The Stellar Mass/Halo Mass relation 



5 




Fig. 4. — The Stellar Mass vs Halo Mass. Black Solid Dots: The SHM relation from our simulations set with stellar masses measured 
using Petrosian magnitudes and halo masses from DM-only runs. This procedure mimics the one followed in M12. Open Dots: Unbiased 
stellar masses measured directly from the simulations. Solid Line: Observational results from M12. 



2.3. The Baryon Content of the Simulated Galaxies 

The more massive galaxies in our sample are disk dom- 
inated, transitioning to irregular galaxies below ~10 10 
M©. The outflows (and the UV cosmic field in halos 
below ~10 9 Mq) significantly lower the baryon fraction 
of the host halos, with a strong trend of lower baryon 
fractions at smaller halo (stellar) mass (see Figure 1). 
However, because the fraction of remaining gas turned 
into stars is low at all galaxy masses and especially in 
dwarfs, the galaxies in our sample have relatively high 
cold gas to stellar mass ratios, typical of real galax- 
ies over the tested range. In Figure 3 we compare the 
cold gas (HI) content of the simulated sample to the 
nearby HI surveys ALFALFA jGiovanelli et al.ll2005h and 
SHIELD (jCannon et al.ll2011l) . The observed dataset in- 
cludes only galaxies closer than the Virgo cluster, for 
a better comparison with our sample of field galaxies. 
With the caveat that selection effects can still play a 
role, there is very good agreement between simulations 
and observations. 

The HI masses are directly measured from the simu- 
lations, where HI and HJ3 abundances are calculated on 
the fly. Magnitudes are measured in the SDSS r-band. 

We verified that the cold gas fraction in the com- 
parison runs adopting a low density threshold for SF 
is about a factor of ten lower at all halo masses. 
Lower resolution simulations of small mass systems 
have often re ported the formation of relatively gas 
poor galaxies (IGoyernato et all I2007t ICoh'n et all 120101: 
lAvila- Reese et al.l 120111 : ISawala et al.l 1201 ID . Low gas 
content in simulated low-z galaxies is likely due to a high 
efficie ncy of gas to stars conversion (|Piontek fc Steinmetzl 
[2011 and/or to an excessive loading factor of the SN 
winds. 

9 note that H2 masses are small and, while neglected in Figure 
3, contribute little to the overall cold gas mass 



3. THE STELLAR MASS - HALO MASS RELATION 

Once individual galaxies in our sample have been iden- 
tified with AHF, the Stellar Mass - Halo Mass ratio 
for each galaxy can then be obtained. The definition 
of "halo mass" includes all DM and baryons within an 
overdensity of 200, but not the mass associated with in- 
dividual satellites (a few % of total at the most). All 
stars not in satellites, but within R200 are associated 
with the central galaxy in the halo. This simple ap- 
proach is si milar to what has b e en done in several pre- 
vious works (jSawala et al.l 120111: IBrook et al.l 120121) and 
similar to what has been used in previous comparisons 
between simul ations and the SHM relation obtain ed us- 
ing the AMT (|Guo et al.H20Tot iMoster et alJl2ljl2l h Our 
new sample of simulated field galaxies (open circles in 
Figure 4) follows closely the shape and normalization 
of the present day SHM relation presented in M12 over 
the 10 9 -10 12 MQ halo mass range (solid line), with 
M s t ar oc M^ alo . This is a large improvement over most 
publis hed works and confirms results on smaller sam- 
ples (IGovernato et al.ll201d:lQh et aHl2011tlGuedes et al.l 
120111: IBrook et al.ll2012h that adopted similar SF and 
feedback recipes. 

Clearly, an approach where a) SF is limited to dense, 
H 2 -rich gas clouds (a highly correlated situation, see 
iChristensen et al.ll2012|) and b) feedback is hydrodynam- 
ically coupled to outflows significantly reduces the SF ef- 
ficiency and the present day stellar mass in galaxy sized 
halos over a wide mass range. Our simulations show that 
both a) and b) are alone not sufficient, but the combi- 
nation is sufficient. We have first verified that in the 
absence of feedback the SF efficiency remains high even 
if the consumption rate of gas is long, as over the course 
of a Hubble time most cold gas within the galaxy even- 
tually turns into stars. Moreover, lack of SN feedback 
fails to remove the low angular momentum gas, origi- 
nating galaxies with an excessive spheroidal component 
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(jGovernato et al.ll20ll iBrook et al.|[20Tll) . Similarly, if 
the identical SN feedback recipe used in this work is ap- 
plied to simulations where SF is allowed in cold, but rel- 
atively low density gas (e.g 0.1 amu, as often adopted 
in the past), it fails to significantly lower the overall SF 
efficiency. In our sample, the mass in stars formed by 
z = in galaxies with H^/high density regulated SF is 
lower by as much as fifteen compared to simulations of 
the same halos adopting a lower density threshold. The 
overproduction of stars in the low threshold runs occurs 
even when metal lines cooling is neglected. 

In summary, while the SF efficiency in the high thresh- 
old simulations is lower, the cold gas content is similar 
to that observed in real galaxies (see §2.3). Hence, a low 
SF efficiency was not obtained by simply increasing the 
feedback strength and "blowing away" all the baryons. 
Combined, these results show that adopting a more re- 
alistic description of where stars form and how feedback 
regulates SF leads to realistic simulations of galaxies. 

However, for a meaningful comparison with observa- 
tions and the SHM inferred in M12, it is important to 
infer both the stellar and halo mass from the simulations 
using the same techniques as the observations. This ad- 
ditional step is necessary as simulations directly measure 
the mass distribution, while observations infer the stellar 
mass from the light distribution. Below we show that this 
more accurate approach affects the results substantially. 
We will provide, for the first time, an accurate compar- 
ison with the present day SHM relation obtained from 
observational data. We used the following procedure: 

• Magnitudes based on the age and metallicity of 
each star particle were derived using the Star- 
burst99 stellar popul ation synthesis models of 
ILeitherer et al.l ([1999) and IVazquez fc Leithererl 
(120051), adopting a Kroupa (2001) IMF. 

• For each simulated galaxy, Pe t rosian aperture mag- 
nitudes (jBlanton et al.l 120011: lYasuda et all 120011 1 
were obtained in the r band. This step is nec- 
essary as observations are limited by the surface 
brightness of the target galaxy dropping below the 
sky brightness. This systematic bias underesti- 
mates the amount of light associated with indi- 
vidual galaxies, and it is estima ted to be of the 
order of 20% for r eal galaxies (|Dalcantonl 119981 : 
Blan ton et al.l 1200 1|) . As our galaxies have light 
profiles that closely mimic those of real galaxies 
(Bro oks et al.ll201ll ). applying this constraint is ap- 
propriate. We verified that the amount of light lost 
is similar to that estimated for observational sam- 
ples. 

• The stellar mass of each galaxy was then estimated 
based on its B-V color and V total magnitude, as- 
suming a Salpctcr IMF to be consis tent with adopt- 
ing the same fitting formula as in IBell fc de Jonsj 
pOOll ). namely L v = ICT^" 4 - 8 )/ 2 - 5 and then 

M star =L V X 10-0.734+1.404x(B-y)_ Wg then ut ;_ 

lize a conversion from Salpeter to Chabrier IMF 
to remain consistent with M12. We find that 
this procedure systematically underestimates the 
true stellar masses (by summing all star particles 
within R V i r not in satellites) of galaxies by 20-30%. 



This result extends over the whole range of galaxy 
m asses. We find t h at the specific criteria adopted 
in IBell &; de Jond (|2001D tends to underestimate 
the contribution from old (i.e., high M/L) stellar 
populations. 

• The halo mass for each simulated galaxy was 
measured re-running each simulation as DM-only, 
matching halos between the two runs and count- 
ing all mass with K V i r (again defined as the radius 
within which the average overdensity is < p >= 
200p cr i t ). This st ep is necessary to fo llow the pro- 
cedure adopted in lMoster et al.l ([2012D . where halo 
masses were obtained from a large cosmological 
simulation that did not include gas physics. This 
procedure avoids a subtle, but significant and sys- 
tematic bias between the total halo mass measured 
in DM-only runs vs those of the same halos in sim- 
ulations that include gas physics and feedback. In 
the latter simulations, feedback can remove a sig- 
nificant fraction of the baryons from the halo, de- 
creasing the total mass within a fixed physical ra- 
dius. The virial radius, if defined at a fixed over- 
density, then shrinks, leading to a smaller total halo 
mass. The decrease in M^j r varies with mass, as it 
depends on the amount of baryons lost, but it can 
be significant, up to 30%. Since the lowest mass 
galaxies in our sample have lost the most baryons 
(in winds and UV background heating), they can 
experience a decrease of ~30% in halo mass. At the 
high mass end, where galaxies retain most of their 
baryons, the simulated galaxies still see a change 
of 5-10% in total halo mass compared to the DM- 
only run (where obviously no baryon mass loss is 
possible) . These results are consistent with esti- 
mates in ()Sawala et al.H2012t) . ISawala et all ()2012[ ) 
also interprets this shift as a systematic reduction 
in the matter infall rate. As even small amounts of 
baryons are removed, the gravitational attraction 
on surrounding material decreases, leading to a de- 
crease in the infall rate of both gas and DM and, 
overtime, to a smaller halo mass. Neglecting this 
effect results in moving the simulations datapoints 
to the left, away from the SHM relation inferred 
using DM-only runs. This bias is particularly no- 
ticeable at small galaxy masses, where the SHM 
relation is steeper. 

In Figure 4, the black solid dots show results from 
our simulations dataset, but using the procedure out- 
lined above, which closely matches that adopted in M12. 
The normalization of the SHM is ~ 40% lower than that 
inferred using the simulations quantities at face value 
and closer overall to the SHM relation of M12. There 
is very good agreement between the SHM inferred from 
the 'artificial observations' and the fit in M12 over at 
least 4 orders of magnitude in halo mass. The mass 
of the brightest galaxies (close in mass and morphology 
to Milky Way analogues) goes from being higher to be- 
ing slightly lower than the SHM. This result confirms 
that the "blastwave" feedback implementation is able to 
reduce the SF efficiency not only at small halo masses 
but also in present day halos around 10 12 M© (see also 
IGuedes et afll20TTI) . 
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Fig. 5. — Top Panel: Halo mass ratio of galaxies in runs with 
baryons and SF vs DM- only runs. Individual halos in DM-only 
runs are typically 30% more massive than their counterparts in 
simulations with gas physics and SF. The effect is smaller in more 
massive halos, whe re baryon loss due to feedback is less (see also 
ISawala ct al. (2012j)). The dashed horizontal line marks the ratio if 
halos had a 100% baryon loss. Bottom panel: Estimated vs. True 
Stellar Mass as a function of halo mass. The stellar mass using ar- 
tificial Pctr osian magnitudes and measured using the photometric 
method in (Bell & dc Jong 2001) vs the "true" Stellar mass mea- 
sured directly from the simul ations. Stellar masses measured using 
the photometric method in (Bell & dc Jong 2001) in combination 
with the flux loss from applying the petrosian magnitudes arc un- 
derestimated by about 50% across the range of galaxy masses in 
our study. 

In Figure 5a we show the total halo mass ratio be- 
tween the simulations that include baryons and SF vs 
the DM-only ones. As discussed above, halos in DM- 
only runs are consistently (and significantly, about 30%) 
more massive. In Figure 5b we show the ratio between 
the stellar masses obtai ned using a combinatio n of Pct- 
rosian Magnitudes with iBell fc de JonsJ (|2001[ ) M/L ra- 
tios (closely following M12) and stellar masses derived 
directly from the simulations. A systematic bias of about 
50% is evident across the whole mass range. The results 
from this section highlight the importance of a careful 
comparison between simulations and observations. This 
approach further reinforces our findings that the model- 
ing of SN feedback greatly reduces the tension between 
the present time SHM relation inferred from the Abun- 
dance Matching Technique for halos with total mass < 
10 12 M© and the predictions from hydrodynamical sim- 
ulations in a cosmological context. Without fine tuning, 
this better agreement also preserves the morphology of 



the galaxies formed, with disk-dominated massive galax- 
ies and low mass irregulars that are gas rich. In future 
work, we plan to extend this approach to higher redshifts, 
using the appropriate criteria to measure stellar masses 
in hig h-z galaxies (e.g.. lPforr et al.ll2012HMaraston et al.l 
I20I3) . 

4. CONCLUSIONS 

We have measured the SHM (stellar mass - halo mass) 
relation for a set of field galaxies simulated in a ACDM 
cosmology and compared it with the redshift zero predic- 
tions based on data from the SDSS and the Abundance 
Matching Technique described in M12. The comparison 
revealed very good agreement in normalization and shape 
over five orders of magnitude in stellar mass. The new 
simulations include an explicit description of metal lines 
cooling and H 2 regulated SF, and SN driven outflows. 
The combination of SF driven by the local efficiency of 
H2 and outflows reduce the overall SF efficiency over the 
whole Hubble time by almost an order of magnitude com- 
pared to older simulations, with resulting M star oc M^ j . 
While a large fraction of baryons is expelled, especially in 
halos smaller than 10 11 M 0) the resulting galaxies have 
an HI content comparable to those inferred by local sur- 
veys, namely ALFALFA and SHIELD. The same galaxy 
set has a cored central DM densi ty distribution, similar 
to observations of real galaxies (|Governato et alJl2012t 
iBrooks fe Zolotovl[20T2l ). 

This agreement between simulations and observational 
data is due to two systematic factors: 1) An implemen- 
tation of SF that relates the SF efficiency to the local 
H2 abundance in resolved star forming regions, result- 
ing in localized feedback that significantly lowers the SF 
efficiency and 2) "observing" the simulations to prop- 
erly compare them to observational estimates of the 
SHM relation. This approach involved creating artifi- 
cial photometric light profiles of the simulated galaxies 
and estimating stellar masses based on aperture magni- 
tudes. Importantly, it also requires coupling the stel- 
lar masses to halo masses derived from DM-only simula- 
tions, rather than the baryonic simulations. Our analy- 
sis shows that adopting photometric stellar masses con- 
tributes to a 20-30% systematic reduction in the esti- 
mated stellar masses. Stellar mass estimates based on 
one band photometric magnitudes are likely to underes- 
timate the contribution of old stellar populations (reflect- 
ing the larger contribution to the total flux coming from 
younger stars). This systematic effect is further exacer- 
bated by the use of aperture based magnitudes, adding 
another 20-30% due to neglecting the contribution of low 
surface brightness populations. Finally, a third system- 
atic effect comes from a difference in halo masses in colli- 
sionless (DM-only) simulations vs simulations including 
baryon physics and outflows. Baryon mass loss makes 
halo masses smaller by up to 30% when calculated at the 
same overdensity (200 in our paper and M12). The effect 
of removing these biases is to move the simulation points 
in Figure 4 further lower and to the right, closer to the 
SHM. 

Notwithstanding the improvements described in this 
and other recent works, further adjustment to our nu- 
merical schemes to model SF and feedback processes are 
most likely required, as more observational constraints 
become available and our understanding of SF improves. 
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In future work we plan to extend our analysis of the 
stellar mass - halo mass relation to higher redshifts and 
larger galaxy masses. The results presented in M12 point 
to a possible discrepancy between the shape of the star 
formation history of real galaxies vs the simulated ones. 
Given the difficulty to obtain robust estimates from faint 
and distant galaxy samples, we expect that the approach 
outlined in this work, i.e. creating artificial observations 
to more directly compare simulations with observations, 
will play an important role. 
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